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Temperature distribution on the receiver surface of a solar power tower plant is of great importance. High 
temperature gradients may lead to local hot spots and consequently failure of the receiver. The temper¬ 
ature distribution can be controlled by defining several aiming points on the receiver surface and adjust¬ 
ing the heliostats accordingly. In this paper, a new optimization algorithm which works based on the 
principles of genetic algorithm is developed to find the optimal flux distribution on the receiver surface. 
The objective is to minimize the standard deviation of the flux density distribution by changing the aim¬ 
ing points of individual heliostats. Flux distribution of each heliostat is found by using the HFLCAL model 
m. which is validated against experimental data. The results show that after employing the new algo¬ 
rithm the maximum flux density is reduced by an order of magnitude. The effects of number of aiming 
points and size of the aiming surface on the flux density distribution are investigated in detail. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

In a solar power tower plant, the receiver plays an important 
role of intercepting reflected solar radiation from the heliostat field 
and transferring thermal energy to the heat transfer fluid. The main 
challenge associated with this process is the high temperature gra¬ 
dient at the receiver surface which may lead to local hot spots, and 
consequently, degradation and failure of the receiver [2] . The tem¬ 
perature distribution on the receiver surface depends on the design 
of the receiver, thermophysical properties of the absorber, heat 
transfer fluid, and the heat flux distribution [3]. Distribution of 
the heat flux on the receiver surface is the only factor which is clo¬ 
sely connected with the performance of the heliostat field. There¬ 
fore, it can be controlled by defining several aim points and 
adjusting the heliostats. 

Two aiming techniques were described by Kistler [4]. One tech¬ 
nique, which is called one-dimensional smart aiming, is to focus 
the heliostats along the height of the receiver until the spillage loss 
starts to increase. The heliostats that are closer to the tower are 
usually focused at the top or the bottom of the receiver while those 
which are farther are aimed closer to the center of the receiver sur¬ 
face. The two dimensional smart aiming is similar, except the 
images are distributed in two dimensions. This technique is usually 
recommended for rectangular cavity apertures or flat plates, as 
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using other shapes of the receivers may lead to increase in the 
spillage loss. Although these methods perform well, the technolog¬ 
ical development of thermal receivers and increasing size of the 
power plants demand more sophisticated aiming strategies. 

In a paper presented by Garcia-Martin et al. [5] an automatic 
closed-loop control method was developed to optimize the tem¬ 
perature distribution within a volumetric receiver at PSA power 
plant. The method is based on measuring the temperature at differ¬ 
ent points on the receiver surface and transferring the power from 
one area to another by changing the aiming points of the heliostats 
when the temperature reaches a maximum tolerance value. In 
another study, Salome et al. [3] presented an open loop approach 
to control the flux distribution on the surface of a flat plate recei¬ 
ver. In this method, a set of aiming points are defined and a grid 
is formed on the surface of the receiver. Then, an optimization 
algorithm called “TABU+ specific neighborhood” is used to find 
the best aiming point for each heliostat. The objective is to mini¬ 
mize the flux spread, F max -F min , while keeping the spillage loss 
above a predetermined value. At the first step all heliostats are 
focused on the center of the receiver. At each iteration, one helio¬ 
stat is selected and its aiming point is changed. If the modification 
leads to an improvement in the objective function, it will be saved 
for the next iteration. It was shown in the paper that the spread of 
the flux density is decreased by 43% with an added spillage loss of 
10% using the proposed algorithm. 

Optimization in engineering design has always been a subject of 
interest to engineers. The genetic algorithm (GA), as one of the 
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Nomenclature 



Am 

mirror surface area (m 2 ) 

p h 

total power sent by a heliostat 

COSW 

cosine factor 

Pc 

crossover probability 

D 

distance between the center of the heliostat and the aim 

Pm 

mutation probability 


point (m) 

P 

1 mm 

second mutations probability 

d 

square root of the heliostat area (m) 

cos rec 

incidence cosine of the reflected central ray from the 

F 

flux density ( ] M) 


heliostat on the receiver surface (°) 

f 

1 max 

\ m ) , X 

maximum flux density f 

W s 

P 

image dimension in sagittal plane 
reflectivity factor 

F 

1 mm 

minimum flux density 

7 HF 

effective deviation (mrad) 

f 

focal distance (m) ; 

7sun 

sunshape error (mrad) 

fat 

atmospheric attenuation factor 

7bq 

beam quality error (mrad) 

H c 

image dimension in tangential plane 

7 ast 

astigmatic error (mrad) 

Id 

direct normal irradiation ) 

7 t 

tracking error (mrad) 

m 

\ m / 

number of heliostats 

7 S 

mirror slope error (mrad) 

n 

number of aiming points 




most popular optimization techniques, has been found very useful 
in solving complex real-world design optimization problems since 
it works with a population of candidate solutions, not a single 
point in the search space. This helps to avoid being trapped in local 
optima as long as the diversity of the population is well preserved 
[ 6 ]. 

Over the past decade, genetic algorithm has been extensively 
used for the optimization of solar thermal systems [7]. Varun and 
Siddhartha used GA to optimize the thermal performance of a flat 
plate solar air heater [8 . Loomans and Visser applied GA for the 
optimization of a solar water heater system [ 9 ]. Godarzi et al. 
[ 10 ] employed GA to optimize the performance of a solar absorp¬ 
tion chiller. GA and artificial neural network were used together 
by Kalogirou 11 to find the optimum combination of the collector 
area and storage-tank size for a solar industrial process heat sys¬ 
tem. Baghernejad and Yaghoubi [12] conducted an exergoeconom- 
ic analysis and optimization of a 400 MW integrated solar 
combined cycle system using GA. Ahmadi et al. [ 13 ] investigated 
multi-objective optimization of a solar dish Stirling engine using 
GA by considering three objective functions, i.e. output power, 
overall thermal efficiency, and rate of entropy generation. Cabello 
et al. [ 14 ] developed a program based on GA to find the optimal 
size of the solar collector area, thermal storage and power of the 
auxiliary system in a direct steam generation power plant. 

In this paper, a new optimization approach based on the princi¬ 
ples of GA is proposed to find the optimal flux distribution on the 
receiver surface of a solar power tower plant. The objective is to 
minimize the standard deviation of the flux density distribution 
on the receiver surface by changing the aiming points of individual 
heliostats. The HFLCAL method [1 ] is used to find the flux distribu¬ 
tion of individual heliostats and is validated against experimental 
data. After presenting the optimization methodology, the final 
results are provided and the influences of different parameters 
are investigated. 


2. Flux density model 

The flux density on the surface of a receiver can be found 
numerically or analytically. In the numerical approach, called 
Monte Carlo ray tracing method, a large number of rays are gener¬ 
ated and traced through different optical stages. A flux map on the 
receiver surface can be generated from the intersection of the 
reflected rays and the surface. SolTrace software [ 15 ], developed 
by the National Renewable Energy Laboratory (NREL), can be used 
to predict the flux density distribution on the receiver surface 


accurately using the Monte Carlo ray tracing method. Although 
the numerical method is very accurate, it requires large computa¬ 
tional time. 

Two well-known analytic models that are used to evaluate the 
flux density and interception efficiency are the UNIZAR model from 
the Universidad de Zaragoza [16] and the HFLCAL model from DLR 
(the German Aerospace Center) [ 1 . According to [ 17 , both of these 
models are appropriate tools though HFLCAL is simpler and slightly 
more accurate than UNIZAR. The HFLCAL model is used in this 
study to evaluate the flux density distribution of each heliostat 
on the receiver surface. The flux map resulting from the entire field 
is generated by superimposing the flux densities of the individual 
heliostats. 


2.2. HFLCAL Model 


In the HFLCAL model a circular normal distribution is used to 
find the flux density distribution on the receiver surface, given as 

[ 3 ]: 


F(x,y) 


Ph 

271(7 HF 



(x - Xf) 2 + (y - y t ) 

2 Gup 




where P h is the total power reflected by a heliostat, g hf is the effec¬ 
tive deviation, and ( x t , y t ) are the coordinates of the aiming point on 
the receiver surface. P h is given as: 

P h =bA m -coswf at -p (2) 

where I D is the direct normal irradiation (™), A m is the mirror area, 
cosw is the cosine factor of the angle between the sun ray and the 
normal to the heliostat surface, f at is the atmospheric attenuation 
factor, and p is the reflectivity of the heliostat. 

The effective deviation, g hf , is the result of the convolution of 
the four Gaussian error functions, namely, the sun shape error 
due to the non-uniform distribution of the solar intensity across 
the sun disk (cr sun ), the beam quality error due to the mirror slope 
error (G bq ), the astigmatic error representing any extra deforma¬ 
tion of the reflected ray if the incident ray is not parallel to the mir¬ 
ror’s normal (cr ast ), and the tracking error (cr t ). g hf is given as [17]: 


'D 2 al 


sun 


+ G bq + G lst 


(7hF 


+ a 2 ) 


y/cos rec 


( 3 ) 


where D is the distance between the center of the heliostat and the 
aim point, and cos rec represents the cosine of the angle between 
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Fig. 1. Contours of flux densities for a heliostat in the PSA power plant: (left) measured and calculated values using the UNIZAR model [17]. (Right) calculated values using the 
HFLCAL model. 
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Fig. 2. A sample form of a candidate solution. 


the reflected ray and the normal to the receiver surface. The beam 
quality error is due to imperfections of the heliostat surface and is 
related to the slope error by: 

(T bci = (2<7 s ) 2 (4) 

The standard deviation of the astigmatic error is given as: 



where H t and W s are the image dimensions in the tangential and 
sagittal plane [18] and are given as: 



d 


D 

—p .— cos w 




D 

7 


cos w - 1 



In the above equations,/represents the focal distance and d is a 
general dimension of the heliostat. In this study, d is equal to the 
square root of the heliostat area. 

The interception efficiency, defined as the fraction of the 
reflected power that reaches the receiver surface at a certain point 
in time, is calculated by [1]: 


dint = 


1 


2nc7 2 HF j x jy 



exp 


(x - x t f + (y -y t ) 2 \ 






dxdy 



A model is developed in MATLAB based on the HFLCAL method 
to calculate the flux density at the receiver of a solar power tower 
plant. The model is validated against the experimental data avail¬ 
able from the PSA power plant [17]. Each heliostat used in the 
experimental study was 6.6778 m wide and 6.819 m high and con¬ 
tained 12 spherical facets (1.105 x 3.010 m). However, all the 


heliostats were modeled as single mirrors in the simulation. The 
total mirror area of each heliostat was 39.9126 m 2 . The informa¬ 
tion about each heliostat including the coordinates, actual distance 
to the receiver plane, and the focal length are given in [17]. It is 
assumed that the direct normal irradiation is 1 kW/m 2 , the mirrors 
are perfect reflectors, and atmospheric attenuation is negligible. 
Therefore, the total power reflected by each heliostat using Eq. 
(2) is given as: 

P h = 39.9126 x cos w (kW) (9) 

Using Eq. (1) and the information given in [17], the flux densi¬ 
ties for the heliostats are obtained. Fig. 1 depicts the contours of 
flux densities for one of the heliostats. The left figure shows the 
measured as well as the simulated values using the UNIZAR model. 
The right figure depicts the flux density distribution obtained from 
the HFLCAL method. As can be clearly seen, analytical methods are 
not able to exactly predict the real shapes of the flux density con¬ 
tours due to the circular symmetry assumption, however, they are 
able to predict accurately the flux density levels on the receiver 
surface. By comparing the two figures, it can be concluded that 
the HFLCAL method works very well and can be confidently used. 


3. Optimization algorithm 

An optimization algorithm is developed based on the principles 
of GA. GA belongs to the larger class of evolutionary algorithms 
that mimic the process of natural evolution in order to find the best 
solution to an optimization problem 19]. The genetic operators 
such as mutation and crossover are applied to a population of 
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Fig. 3. Crossover operation on two of the candidate solutions. 
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as shown in Fig. 2. The heliostats are numbered from 1 to m. Sim¬ 
ilarly, the aiming points are numbered from 1 to n. The matrix 
shown in Fig. 2 indicates that the first heliostat is focused on the 
first aiming point; the second heliostat is focused on the third aim¬ 
ing point; and the mth heliostat is focused on the (n - 1 )th aiming 
point. 

The number of matrices that are randomly generated in this 
step, called the number of initial population, is defined by the user. 
A MATLAB cell data structure is used for the programming [20]. 


Fig. 4. Mutation operation on one of the candidate solutions. 


3.2. Crossover operator 



X(m) 


Fig. 5. Flux density map for a randomly generated binary matrix. 

candidates. The best candidates are selected at the end of each gen¬ 
eration by measuring the fitness function. The process is repeated 
for a number of generations until the best solution is found. The 
optimization algorithm can be summarized as: 

(1) Initial population is randomly generated. 

(2) The fitness function is evaluated for each population. 

(3) The offsprings are generated by applying genetic operators, 
i.e. crossover and mutation, to probabilistically selected 
individuals from the initial population. 

(4) The fitness function is evaluated for the offsprings. 

(5) New population is selected and the algorithm continues 
from step 2. 

In this paper, the genetic algorithm is employed to find an aim¬ 
ing strategy for a solar power tower plant in order to distribute the 
flux uniformly on the receiver surface. The main objective is to 
minimize the standard deviation of the flux density distribution, 
which is found by measuring the flux density at multiple points 
on the receiver surface. The detailed methodology is explained in 
subsequent sections. 

3.1. Generating initial population 

In the GA, each candidate is usually represented as a vector con¬ 
taining information about all the design variables. In this problem, 
each candidate is represented by a matrix where the numbers of 
columns and rows are equal to the numbers of heliostats and aim¬ 
ing points, respectively. Therefore, for n number of aiming points 
and m number of heliostats, the size of the matrix will be n x m, 


In the crossover process two of the candidate solutions are 
selected and combined to produce children, which are the new 
candidate solutions. The crossover probability (P c ) is the parameter 
that controls how often the crossover will be performed, which is 
defined by the user [19 . The crossover process consists of three 
main steps as shown in Fig. 3: 

(1) Two of the candidate solutions are randomly selected. 

(2) A column number is randomly selected. 

(3) The values are swapped between the two matrices following 
the selected column number. 

The main intention of using crossover is to find new solutions 
that contain good parts of the old solutions; however, it does not 
necessarily mean that the offsprings are any better than the 
parents. 

3.3. Mutation operator 

Mutation is another genetic operator which helps to avoid being 
trapped in the local minimum and maintains the diversity in the 
population [20 . Similar to crossover probability, mutation proba¬ 
bility (P m ) controls how often the mutation is performed on the 
old population. The mutation is incorporated into the program by 
the following steps: 

(1) A candidate solution from the old population is randomly 
selected. 

(2) One of the columns of the chosen matrix is selected. The ele¬ 
ment in that column that is equal to 1 is changed to 0. 

(3) Another element of that column is randomly selected and is 
changed from 0 to 1. 
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Fig. 6. The receiver surface, aiming surface, and the generated squares that are used 
to explain the algorithm. 
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Considering the physics of the problem, it is expected that per¬ 
forming mutation only on one column will not change the objec¬ 
tive function very much. In other words, changing the aiming 
point of a single heliostat will not affect the overall standard devi¬ 
ation of the flux density distribution significantly. Therefore, the 
program is written in a way that allows mutation to be performed 
on multiple columns of the selected matrix. The number of muta¬ 
tions is defined by the second mutation probability (P mm ). In other 
words, for each column a random number is generated and muta¬ 
tion will be performed on that column if the random number is less 
than P mm . 

Fig. 4 depicts the mutation operation on a randomly selected 
matrix. As it can be seen, mutation is performed on the second 
and (m - l)th columns, which means that the aiming points of 
the second and (m - l)th heliostat are changed. 

3.4. Selection 

The last step is to select the best candidates from the old popu¬ 
lation and the offsprings for the next generation. Tournament 
selection strategy is used for this purpose which involves running 
several tournaments between two randomly chosen candidates. 


The winner of each tournament is the one with the lowest standard 
deviation of the flux density distribution on the receiver surface. 
The tournament competitions continue until the number of the 
selected candidates equals the number of in initial population. 
One generation is completed by selecting the new population 
and the program proceeds to the next generation. The program 
continues to run until the maximum number of generations 
defined by the user is reached or some other termination criterion 
is met. 

3.5. Self-modifying algorithm 

As it is already discussed, the optimization algorithm is based 
on generating random binary matrices. The flux density maps pro¬ 
duced by these random matrices do not usually reflect the 
expected shapes. In other words, the maximum flux density might 
be located at the corner of the receiver or the contours may not fol¬ 
low a reasonable pattern. This problem is clearly shown in Fig. 5. 
Therefore, another algorithm, called self-modifying algorithm, 
needs to be developed to modify the flux density maps and turn 
them into more realistic shapes by changing the elements of the 
random matrices. In this program the ideal map is defined as the 
one in which the maximum flux density is located as close as 
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Fig. 7. Flowchart for the self-modifying algorithm. 


























S.M. Besarad et al.J Energy Conversion and Management 84 (2014) 234-243 


239 



Fig. 8. Modification of flux density map shown in Fig. 5 by changing the aiming points of the heliostats (a) 3rd iteration, (b) 8th iteration, (c)16th iteration, (d) 24th and the 
last iteration. 


possible to the center of the receiver and gradually decreases as it 
gets closer to the receiver sides. 

The algorithm developed to handle this task can be explained 
with the aid of Fig. 6. The largest square represents the receiver 
surface and 81 aiming points are distributed on a smaller square, 
which is called aiming surface. 

The discussion about the size of this surface and the number of 
aiming points is presented later. As shown in Fig. 6, the aiming 
points can be connected to each other to form new squares. Begin¬ 
ning from the squares #1 and #2 as the outer and inner squares, 
respectively, the flux densities for all the aiming points located 
on the sides are calculated. The maximum flux on the square #1 
should be less than the minimum flux on the square #2. If this con¬ 
straint is violated, one of the heliostats that are focused on the aim¬ 
ing point representing the maximum of square #1 is repositioned 
and is focused on the point having the minimum flux density of 
square #2. This algorithm continues until the flux density at all 
the points on square #1 becomes less than those on square #2. 
Next, the outer square is replaced by square #2 and the inner 
square is represented by square #3 and the process continues. 
The flowchart for the algorithm is shown in Fig. 7. 

During optimization, all the generated matrices are self-modified 
by this algorithm. The algorithm is applied to the flux map shown in 
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Fig. 9. Heliostat field for a 50 MWth solar power tower plant in Daggett, California. 
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Table 1 

Field parameters. 

Heliostats 


Width 

12.84 m 

Height 

9.45 m 

Reflectivity 

0.88 

Receiver 


Tower height 

115 m 

Tilt angle of the aperture 

12.5° 

Aperture width 

12 m 

Aperture height 

12 m 




5000 kW 
m 2 
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Fig. 10. Flux density map for the solar power tower plant assuming all the 
heliostats aim at the center of the receiver. 
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Fig. 11. Distribution of the aiming points on the receiver surface. 


Fig. 5. Fig. 8 depicts how the flux density map is modified with the 
number of iterations. As can be seen, on the 24th iteration the max¬ 
imum flux density is located at the center and gradually decreases 
until it reaches the minimum level close to the sides of the receiver. 

4. Results and discussion 

4.1. Optimization results 

As a case study, a field layout as shown in Fig. 9 was designed 
for a 50 MWth system in Daggett, California. There are 580 helio¬ 



stats which are placed on the field based on the biomimetic layout 
proposed by Noone et al. 21 . This field and the methods explained 
in the preceding sections are employed to find the optimal flux 
density distribution on the receiver surface. Other field parameters 
are given in Table 1. It is assumed that a sun is 2.51 x 10~ 3 mrad, o s 
is 10 -3 mrad, the tracking error is negligible, and the facet canting 
is on-axis parabolic. 

On March 21st at noon, the direct normal radiation in Daggett is 
given as 0.926 ^ [22]. Assuming all the heliostats aim at the cen¬ 
ter of the receiver at this time of the year, the flux density map 
using the FIFLCAL method is shown in Fig. 10. As can be seen, the 
flux density at the center of the receiver is close to 5000 which 
is extremely high. Therefore, an aiming strategy is required to uni¬ 
formly distribute the flux over the entire surface. 

In this regard, a number of aiming points are distributed on the 
receiver surface, as shown in Fig. 11. All of these points are located 
inside a square which is called the aiming surface. Then, the opti¬ 
mization algorithm is used to find the appropriate aiming point for 
each heliostat that leads to the least standard deviation of the flux 
density distribution on the entire surface. In this study, the aiming 
surface is a 8 x 8 square on which 81 aiming points are uniformly 
distributed. The evolutionary process of optimization is accom¬ 
plished with a population size of 100 with crossover probability, 
P c , mutation probability, P m , and second mutation probability, 
P mm , of 0.9, 0.2, and 0.5 respectively. Fig. 12 shows the variations 
of the objective function with the generation number. The optimi¬ 
zation algorithm is terminated after 235 generations when the 
objective function reaches a plateau. 

Fig. 13 depicts the optimal flux distribution on the receiver sur¬ 
face. It can be seen that the maximum flux density is 550 kW/m 2 , 
which is almost one tenth of the maximum flux when all the helio¬ 
stats aim at the center of the receiver (Fig. 10). Therefore, it can be 
concluded that the new algorithm is able to successfully distribute 
the flux on the entire surface. The interception efficiency using the 
proposed aiming strategy is 0.9214 as opposed to 0.9906 in the ori¬ 
ginal case. According to Fig. 13, there is a spot which is shown as 
dark red and is not located at the center of the surface. The reason 
is that the self-modifying algorithm is based on calculating and 
comparing the flux densities only at the aiming points. Therefore, 
there might be some points in between that cannot be captured 
by the algorithm. This is not a concern as long as the size of the 
spot is small. Flowever, increasing the number of the aiming points 
can fix this problem, which is discussed in the following sections. 
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A map of the heliostat field is presented in Fig. 14 to show 
which heliostats were aimed at which points. The numbering for 
the aiming points begins from the bottom left corner in Fig. 11 
and ends at the top right. 


4.2. Influence of the aiming surface size 


The size of the aiming surface is an important parameter which 
directly affects the flux density distribution as well as the intercep¬ 
tion efficiency. A smaller aiming surface leads to higher intercep¬ 
tion efficiency and higher maximum flux density, and vice versa. 
Fig. 15 shows the variations of the interception efficiency and the 
maximum flux density with respect to the size of the aiming sur¬ 
face. In this case, the aiming surface is represented by a square 
on which 81 aiming points are distributed. The calculations were 
performed for March 21st for Daggett, California. 

As can be clearly seen in Fig. 15, the smallest size of the aiming 
surface (5 m x 5 m) has the highest interception efficiency and the 
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Fig. 16. Optimal flux density distribution on the receiver surface: (left) 49 aiming points (right) 121 aiming points. 


maximum flux density increases to 1070 kW/m 2 . As the size of the 
aiming surface increases to 12mx 12 m, the interception effi¬ 
ciency decreases to 0.7533. This means that close to 25% of the 
energy reflected by the heliostats is lost and not intercepted by 
the receiver. This has a huge impact on the performance of the 
power plant. 

Therefore, one has to consider the maximum heat flux density 
that can be tolerated by the receiver surface and size the aiming 
surface accordingly. It is noteworthy that by decreasing the size 
of the aiming surface, it might be necessary to decrease the num¬ 
ber of aiming points as distributing a large number of aiming 
points on a small surface can be a challenge from practical point 
of view. 

4.3. Influence of the number of aiming points 

Increasing the number of aiming points leads to better 
distribution of the heat flux on the receiver surface. Moreover, 
the self-modifying algorithm is more efficient when a large num¬ 
ber of aiming points is employed. However, in a real power plant 
the number of aiming points that can be distributed on a surface 
depends on the size of the surface and uncertainty in aiming the 
heliostats. Therefore, having a very large number of aiming points 
on a small receiver may not be practical nor necessary. 

Fig. 16 depicts the flux maps on the receiver surface when 49 
and 121 aiming points are considered. All other parameters are 
same as those used in Fig. 13. As can be clearly seen from the left 
plot, when there are 49 aiming points, the maximum flux density is 
located on the left corner of the square. However, the size of this 
spot is small and the flux density on this point is not very different 
from its neighbors. Moreover, the flux density is well-distributed 
over the entire surface of the receiver. Therefore, it can be con¬ 
cluded that although the number of aiming points has decreased 
from 81 to 49, the algorithm is still able to distribute the flux very 
efficiently. The same analysis is carried out for 121 aiming points 
and the resulting flux map is shown on the right plot. According 
to the figure, there is no dark spot and the flux density is perfectly 
distributed on the surface and the maximum flux density is about 
500 ^ 

DUU m2- 

From the results shown in this section, it can be concluded that 
it is desirable to have a larger number of aiming points from com¬ 
putational point view, even though the program is still able to dis¬ 
tribute the flux efficiently using a smaller number of aiming points. 
The selection of the number of aiming points depends on the size 
of the receiver and uncertainties in aiming the heliostat. 


5. Conclusions 

The work in this paper presents an optimization algorithm that 
can be used to obtain uniform heat flux density distribution on the 
receiver surface of a solar power tower plant. The flux density of 
individual heliostats was modeled using the HFLCAL method, 
which was validated against experimental data. The optimization 
code was developed based on the principles of GA and modified 
by considering the physics of the problem. The results showed that 
using the new algorithm the maximum flux density is reduced by a 
factor of 10, reaching around 500 The analysis has shown that 
the size of the aiming surface has a huge impact on the intercep¬ 
tion efficiency and the maximum flux density and should be cho¬ 
sen as small as possible. This, of course, depends on the 
maximum flux density that the receiver material can tolerate. 
Moreover, it is shown that the algorithm works very well for differ¬ 
ent numbers of aiming points, however, it is desirable to increase 
this number as much as possible. It is noteworthy that in this paper 
the heliostats are modeled as single mirrors, but the algorithm 
would also be applicable if the heliostats were represented by dif¬ 
ferent facets. 
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